Extrapolation Method for the No-Core Shell Model 



H. Zhan, A. Nogga, and B. R. Barrett 
Department of Physics, University of Arizona, Tucson, Arizona 85721 

J. P. Vary 

Department of Physics and Astronomy, Iowa State University, Ames, Iowa 50011 



o 
o 

(N 



p. Navratil 

Lawrence Livermore National Laboratory, L-414, P-0. Box 808, Livermore, California 94551 

(Dated: February 8, 2008) 

Nuclear many-body calculations are computationally demanding. An estimate of their accuracy is 
often hampered by the limited amount of computational resources even on present-day supercomput- 
ers. We provide an extrapolation method based on perturbation theory, so that the binding energy 
of a large basis-space calculation can be estimated without diagonalizing the H[amiltonian in this 
space. The extrapolation method is tested for '^H and ''Li nuclei. It will extend our computational 
abilities significantly and allow for reliable error estimates. 
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II. NO-CORE SHELL MODEL 
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In recent years a great deal of progress has been made 
in solving the nuclear many-body problem based on mi- 
croscopic nuclear interactions. Using a stochastic ap- 
proach 0, S 01 or directly diagonalizing the Hamiltonian 
01 , previous investigations have been able to improve the 
accuracy of the calculations to a level, so that rather 
small parts of the Hamiltonian, like the three-nucleon 
(3N) forces, can be probed by a comparison of the pre- 
dicted spectra to the experimental values. For unam- 
biguous conclusions, a reliable error bound for the calcu- 
lations is of highest importance. Currently, error bounds 
have often been established using benchmark calculations 
P^IEIEQ' this paper, we propose a new scheme to 
extend no-core shell model (NCSM) calculations beyond 
their current limits, and to estimate error bounds for ex- 
isting calculations. We establish a correlation between 
the expectation value of the Hamiltonian (H) with re- 
spect to the approximate ground state and that of the 
square of the Hamiltonian (i/^), following the ideas of 
Mizusaki and Imada 0, . This will improve the esti- 
mate of binding energies and provide more reliable error 
estimates. 



To this aim, we proceed with a brief introduction to 
the main concepts of NCSM calculations in Section ^ 
Then we motivate the correlation by the properties of 
NCSM effective interactions and perturbation theory in 
Section ITTll A comparison between the extrapolated re- 
sults and the known binding energies for '^H and ^Li is 
given in Scction llVl Conclusions and outlooks to the ap- 
plications of this method and a discussion of the relation 
of this work to previous ones in Refs. 00 are presented 
m Section El 



Usually shell-model calculations assume an inert nu- 
clear core. Taking only the valence nucleons as active 
particles clearly has the advantage of reducing the num- 
ber of many-nucleon states. But, so far, this also means 
that effective shell model interactions have to be used, 
which cannot simply be related to nucleon-nucleon (NN) 
interactions, as they have been developed for few-nucleon 
systems. 

The NCSM approach is different. All nucleons are 
taken to be active, and the same interactions are used 
as they are in traditional few-body methods. The cal- 
culations are performed in a finite antisymmetrized har- 
monic oscillator (HO) basis, using either Jacobi ^3 or 
Cartesian coordinates like those in usual shell model in- 
vestigations. In general, the short range repulsion of 
nuclear forces cannot be easily described in a finite ba- 
sis. In particular, the HO basis is not well suited to 
describe the short range correlations and, on the other 
hand, the exponential tail of bound-state wave functions. 
Consequently, effective interactions appropriate to the 
basis-size truncation must be derived from the under- 
lying nuclear forces in order to achieve convergence with 
a manageable number of basis states. These effective 
interactions can be systematically related to the "bare" 
NN interactions II [llllllll- The scheme described 
in Refs. 0, Il3l | is based on unitary transformations 
of the Hamiltonian , which decouple the model space 
from the complete Hilbert space describing the quantum 
mechanical system. 

The starting point of all NCSM calculations is a non- 
relativistic A-body Hamiltonian, which includes two- 
body interactions. The extension to 3N forces has been 
introduced in [isl IT^ . We do not take them into ac- 
count in this study. Adding the center-of-mass (CM) 
HO potential, which is subtracted at a later stage in the 
calculation, one can cast the Hamiltonian into the 
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where Pi (r^) is the momentum (position) of the particle 
i with mass m^, and = the total mass of 

the A particles. 

It is possible to establish a unitary transformation 
of the Hamiltonian, which decouples two parts of the 
Hilbert space - a rather small finite model space P and 
the rest of the Hilbert space Q. The projection operators 
on the two spaces are also called P and Q. They fulfill 
the relations Q = 1 - P, and Qe~^ H^e^ P = 0, where 
is the unitary transformation 

The solution of the A-body problem is made possible 
by the observation that, for nuclear problems, a cluster 
approximation to the full unitary transformation suffi- 
ciently speeds up the convergence compared to the bare 
interactions, so that practical calculations are possible. 
Instead of solving the full A-body problem, the proce- 
dure is carried out for a much smaller a-body problem, 
with a = 2 or a = 3 in practice. For such cases, one 
starts with the Hamiltonian 
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Note that the mass of the full yl-body system enters 
into the strength constant of the relative HO potential. 
Therefore, for this a-body Hamiltonian, the HO interac- 
tion does not cancel in the relative motion, and provides 
a confining mean field interaction. However, the HO in- 
teraction is cancelled in the A-body calculation. This 
procedure improves the convergence of our A-body re- 
sults with increasing model space size. 

Using the unitary transformation satisfying 
QaB'^'^H^e^'^Pa — 0, one determines an effective 
Hamiltonian H^ffa = PaS'^^H^^^^Pa, which exactly 
describes the a-body cluster in the model space. The 
truncation Pa from the full a-body Hilbert space is 
related to P by the requirement that the a-body states 
included in P are also included in Pa- One then defines 
the effective interaction to be used for the A-body 
calculation as 
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With increasing size of the model space, the effective in- 
teraction converges to the "bare" interaction, so that the 
"bare" problem is recovered, meaning that the approxi- 
mation is controllable. 



Due to the cluster approximation, for a < A we no 
longer have an exact effective interaction. This shows up, 
for example, in an dependence of the binding energy. 
However, experience shows that shell-model calculations 
converge much faster, if performed with effective forces, 
as defined above All calculations in this paper are 
based on effective interactions obtained from two-body 
cluster solutions. 



III. CORRELATION BETWEEN Ea AND A£ 

The binding energy Eq — (H) evaluated in an approx- 
imate ground state must approach the exact binding en- 
ergy £o as the energy variance AE^ = {H"^) — (H)'^ van- 
ishes. One can easily extrapolate the exact binding en- 
ergy with a series of approximate calculations, once the 
behavior of Eq as a function of AE is determined. It 
is suggested by Mizusaki and Imada that Eq can be ex- 
anded in terms of AE^Eq'^ for small values of AE^E^'^ 
[. They propose two extrapolation formulae for tra- 
ditional shell model calculations: Eq ~ a AE^Eq'^ + 8o 
and Eq ~ gq AE^ + aiAE* + Eq, where a, ao, and ai 
are fitting parameters. This has motivated us to search 
for the correlation between Eq and AE in context of the 
NCSM, and establish an extrapolation method for the 
NCSM. 



A. Notation 

For the following investigation we need to explicitly de- 
fine the model spaces. To this aim, we truncate the full 
Hilbert space spanned by the antisymmetrized HO basis 
at a maximum total HO quantum number . This 
counts the number of oscillator quanta including the low- 
est oscillator configuration for the nucleus of interest [23|. 
This truncation ensures that all states are included up 
to a given energy, so that spurious CM motion can be 
projected out. The goal of the following study is to ac- 
quire results for a fixed Nm from calculations for even 
smaller subspaces P oi P truncated by Nm < Nm- The 
effective interactions will be obtained for the larger Nm 
in all cases. 

We now switch to matrix notation, because it helps to 
demonstrate the structure of the effective Hamiltonian, 
H, which can be decomposed into 
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The boundaries of the blocks in H are given by the 
truncation Nm of the subspaces P. H acts only in P and 
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H in the remainder space P — P. Because the effective 
interactions and H are defined for the complete P space, 
Nm completely determines H, so one may write explicitly 
H(A'm). We assume that ip and 4> are the eigenvectors of 
H and H, respectively, i.e., Hipi — Eiipi, i — . . .n — 1, 
and H(/)j = Uj<j)j, j = 0.. 
the dimensions of P and P 
eigenvectors of Hq are 



. n — 1 , where h and n are 
- P, respectively. Then the 
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because of the form of Hq. The ground state of Hp, i.e., 
Aq, is of the from [■(/;J,0]'^, and we denote the ground- 
state energy by Eq. 

The energy dispersion AE is defined as 

Ai;2=AjH2Ao-(ATHAo)2=VjBBTVo, (7) 

where we have made use of the specific forms of Hq and 
Hi. The quantity AE'^ measures how well Aq approx- 
imates the ground state of H. As approaches N,n, 
AE will vanish, and Eq will become the true binding 
energy Sq. 

To establish a scaling between Eq and AE'^, we try to 
improve the results of H using perturbation theory. Sev- 
eral more approximations will be necessary to motivate 
the scaling. These approximations will be tested numer- 
ically in model calculations for the '^H binding energy. 
We do not intend to improve on the currently available 
binding energy results J., i^, J£j JH] , but restrict our- 
selves to the P space truncated at = 20, which yields 
realistic, but not fully converged results. Here we will 
only aim to recover P-space results with Nm < 20 calcu- 
lations. Note that this is in contrast to the usual NCSM 
calculations where one investigates the dependence on 
Nm itself, and the effective interactions are renormalized 
for each Nn ■ This also means that our results will explic- 
itly depend on Q. But the small model space sizes used 
will allow us to extract intermediate results, which will 
demonstrate the origin of the scaling behavior. We will 
also show results for ^Li to demonstrate the applicability 
to more complex systems. 



B. Application of Perturbation Theory 

With the block-block decomposition Eq. Q), second- 
order perturbation theory gives the eigenenergy of H, 
i.e., £i, corresponding to Ei, as 
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i = Q...n-l. 



(8) 



The next non-vanishing perturbative term is of the 4th 
order. Despite its simplicity and accuracy (see Fig. 
Eq. (jHJ is not useful in establishing a correlation between 
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FIG. 1: Perturbative calculations of the binding energy Eq 
of the tritqn with = 20 and hO. = 28 MeV. The values 
of Eq for H are denoted by open squares. The circles are 
calculated using Eq. ||HJ, the downward triangles Eq. 1101 . 
and the upward triangles Eq. 



Eq and AE, because it requires the knowledge of all the 
eigenvalues and eigenvectors of H. 

One may avoid diagonalizing H by a redefinition of the 
decomposition: 
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where D is a diagonal matrix containing the diagonal 
elements of H, and C = H - D. For NCSM Hamilto- 
nians, this decomposition is appropriate (see subsection 
IIII ("p . We note that Eq. 10 remains unchanged. Similar 
to Eq. second-order perturbation theory gives 



E, - D, 



i . . . n — 1 , 



(10) 



where Bj is the j-th column of B. Up to the third order, 
we have 



(2) 
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Figure shows the binding energy of •^H calculated 
using equations JHl, pt)|l . and Ullfl . and compares them 
to the results from the diagonalization of H. This is 
done for two different nuclear interactions, CD-Bonn plj 
and AVIS |23|, using a Jacobi basis. This basis is ob- 
tained from the eigenstates of the antisymmetrizer (see 
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(a) Hamiltonian 




FIG. 2: (a) The effective Hamiltonian of tlie triton witli the 
CD-Bonn potential and li^l = 28 MeV. The basis is sorted ac- 
cording to the quantum number A'^, but the ordering within 
each group of states with the same A*" is arbitrary. This or- 
dering will be assumed in all the following figures, (b) The 
effective kinetic energy, which is calculated in the same way 
as the Hamiltonian, but with the NN interaction turned off. 



[13 )■ The antisymmetrized states have a well defined 
total quantum number N. Within each group of equal 
N, the states are ordered arbitrarily. Terms acting on 
the CM do not contribute in this case. It is seen that 
the perturbative calculations can greatly improve small 
Nm results, regardless of the choice of the nuclear in- 
teraction. Figure ^ also indicates that the third-order 
term in Eq. is small for Nm > 12. This suggests 
that the non-diagonal matrix elements of H, i.e., C, are 
small relative to their associated energy denominators in 
second-order perturbation theory. 



C. Eq vs. AE^ Scaling 

The expected scaling behavior can be motivated by 
considering the features and results of equations H1U|I and 
To this aim, the behavior of Djj has to be under- 
stood in more detail. Figure |21 shows the matrix elements 
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FIG. 3: The diagonal elements of the Hamiltonians and ki- 
netic energies for "^H. The horizontal axis in the upper panel 
enumerates the basis states for each individual diagonal ele- 
ment (crosses) of the model with the CD-Bonn potential. The 
horizontal axis in the lower panel indicates the HO quantum 
number N for other data. Bare T (squares) is the expectation 
value of the kinetic energy operator with respect to HO basis 
states. The remaining results are the averages of diagonal el- 
ements with the same HO quantum number N. From right to 
left, the first step of the individual diagonal elements in the 
upper panel corresponds to A'^ = 20, the next A'^ = 18, and so 



of the effective two-body Hamiltonian and the effective 
two-body kinetic energy for our triton model both in the 
3N basis. The effective kinetic energy is calculated in 
the same way as is the effective Hamiltonian, except that 
the NN interaction is turned off. In this case, the CD- 
Bonn potential is used with Nm = 20 and H^l = 28 MeV. 
Comparing the kinetic energy with the full Hamiltonian, 
one sees that the kinetic energy dominates over the con- 
tributions from the NN interaction, especially for large 
HO quantum numbers. In the figure, the basis states are 
ordered according to their total HO quantum number TV. 
The ordering within each of these groups is arbitrary. 

The potential energy affects mostly the low-iV states 
and the cross terms among them and high- lying states. 
The diagonal elements of the Hamiltonian arc dominant 
over off-diagonal elements, which confirms our expec- 
tation. It is also striking that a large number of the 
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FIG. 4: The linear relation between Eq and AE^. The upper 
panel is for '^H with the CD-Bonn potential, the lower one 
for AV18. Symbols from right to left correspond to Nm ~ 
10, 12, . . . , 20 for each model. The lines fit the results from 
iVm = 10 to 18. 



diagonal elements are roughly equal. The dark blocks 
next to the diagonals of the Hamiltonian and the kinetic 
energy are due to the fact that the kinetic energy op- 
erator changes the quantum number N by ±2 (also 0, 
which corresponds to the low-amplitude diagonal blocks). 
Since everything is renormalized to effectively include the 
higher-space (TV > Nm) influence, ± 4, ±6, . . . terms 
show up in the Hamiltonian and the kinetic energy, and 
they become progressively weaker towards low N . The 
above properties persist for other interactions, and other 
values of ft or Nm- They are a common feature of NCSM 
effective Hamiltonians, at least for the 3N system. 

The behavior of the diagonal matrix elements is quan- 
tified in Fig. 121 The upper panel shows individual diag- 
onal elements, and, again, the basis states are ordered 
according their HO quantum number. The lower panel 
plots the averages of the diagonal elements that have 
equal N. The step structure among the individual di- 
agonal elements reflects the fact that Djj is roughly the 
same in each group of states with the same HO quantum 
number. The reason for the flattening for larger N is 
two-fold. First, the number of states increases dramat- 
ically for larger N and, second, the renormalization of 
the effective matrix elements reduces large N diagonal 
elements. This is seen by comparing the average effec- 
tive kinetic energy (circles) with the bare kinetic energy 



TABLE I: The ground state energies, £o, of H for ^H, their 
deviation from the correct result SE, and the computational 
time necessary for the solution. The results of the small space 
(H) solution are compared to the extrapolation method and 
the perturbative calculation using Eq. IIH . AU results are 
for CD-Bonn potential and fin = 28 MeV. 
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''5£; = _Eo - £0, in units of keV. 

■^In units of tg, the time needed to solve the ground state of H 
with Nm = 8. 

''The extrapolation uses Eq and AE^ of H{iV) with TV from 10 to 

Nm- 



(squares), which is the expectation value of the kinetic 
energy operator with respect to HO basis states. For 
N < 14, both follow the linear ^^^^^ hQ behavior as ex- 
pected. For higher N, the effective kinetic energy turns 
flatter. The flgure also compares the averaged diagonal 
elements of the full Hamiltonians (triangles) to the ones 
of the kinetic energy. This again demonstrates the dom- 
inance of the kinetic energy for > 0. One can expect 
that the same behavior holds for more complex nuclei. 

This structure of the effective Hamiltonian guarantees 
that the denominators of the second-order term in equa- 
tions ((Tnt and l(TT)) are, to a high accuracy, equal over a 
wide range of N. Thus, we have, approximately, 

1 1 

£oc^Eo--J2 {^^oBjY = Eo- -^'^BB'^i^o, (12) 

where a is a positive constant. Since the absolute value 
of Eq is much smaller than Djj, the constant a should 
be roughly the average value of the diagonal elements 
with N > Nm- This also means that is only weakly 
dependent on Eq. Neglecting higher-order perturbative 
terms and taking a constant true binding energy £q, it 
follows that Eq oc AE"^. This motivates a linear scaling 
behavior between Eq and AE^. 



IV. APPLICATION OF THE EXTRAPOLATION 

Now that the linear scaling between Eq and AE'^ has 
been motivated, one can estimate the true binding en- 
ergy £q by a linear regression of £^0 and AE'^ which are 
calculated with Nm < Nm, and extrapolating it to the 
point where AE^ = to estimate £q at A^^. 
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TABLE II: Errors in perturbative estimates of the '^H bind- 
ing energy using equations Q and and the difference 
between extrapolated binding energies and corresponding per- 
turbative results from equation llUl . All results are for CD- 
Bonn potential and ftfi — 28 MeV. 
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"Perturbative estimates are calculated with Nm = Nm — 2. The 
results are in units of keV. 

''£o,e is the extrapolated binding energy using only the results of 
Eq and AE^ from Nm = A^m - 2 and Nm - 4. 





1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 


-15 


■_ %[ Af,„ = 12 

/ti i lO iViC V 




-20 






-25 






-30 


1 1 1 1 1 1 1 1 


□ CD-Bonn - 

° AV18 : 



500 1000 1500 



(McV2) 

FIG. 6: The same as Fig.QJ but for the ^Li nucleus. Symbols 
from right to left correspond to Nm = 6, 8, 10, and 12 for each 
model. The lines fit the results from Nm ~ 6, 8, and 10. 



The numerical results for the relation between Eq and 
AE"^ are shown in Fig. ^ for '^H with different NN in- 
teractions and different values of hil. Nm is 20 in all 
cases, and Nm varies between 10 and 20. Additionally, 
linear fits to the results for Nm = 10 to 18 are plotted. 
Clearly, the linear scaling behavior is confirmed by these 
calculations. 

To be more quantitative, we compare results of direct 
calculations, the extrapolation based on the scaling be- 
havior, and perturbative estimates based on Eq. (|ll|l in 
Table 13 The table indicates the computational efforts 
necessary by giving run times for different calculations. 
One sees that a stable extrapolation is possible starting 
from Nm — 14. The extrapolation error for larger cal- 
culations is comparable to the error of the perturbative 
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FIG. 5: The time consumed in the calculation in units of ts, 
the time needed to solve the ground state of H with Nm = 
8. The diamonds mark the time it takes to extrapolate the 
ground state of H using all the ground states of il{N) with 
A'' = 10 . . . Nm, so it is the accumulation of the time needed 
to solve H(A'^). All other symbols are as in Fig. Q 



estimates, indicating that both can be traced back to 
higher order terms in the perturbative expansion. In cal- 
culations for Nm = 30, we have confirmed that the range 
of the linear behavior is extended for larger Nm ■ This is 
expected because ^/Djj is driven by the kinetic energy 
and, therefore, is proportional to 1/N (see Fig.O. Con- 
sequently, if N is increased by one step (i.e. 2 units), then 
the change in 1/Djj is of order 1/N^, which decreases 
with A^. We also note that the constant a = 268 MeV 
for CD-Bonn with hil = 28 MeV is comparable to the 
diagonal elements Djj shown in Fig. O 

The effect of Nm is demonstrated in Table HTl where we 
list errors in perturbative estimates of the binding en- 
ergy using equations © and ifTUIl . The two equations are 
both of the second order, and the errors tend to decrease, 
though not monotonically, as the model space increases. 
At the same time, Table Hll shows that the extrapolated 
binding energies converge to the results of equation ifTUIl . 
This is expected, because our extrapolation method, i.e., 
equation (|12|l . is based on an approximation of second- 
order perturbation theory. The behavior of the pertur- 
bative calculations and extrapolations suggests that one 
can reduce the extrapolation error by increasing Nm. 

The run times are also compared in Fig. [S] It is clear 
that pcrtmbation theory based on Eq. © does not im- 
prove the timings because of the extra time needed to 
diagonalize H. The extrapolation method for Nm = 14 
yields sufficiently accurate results, but reduces the CPU- 
time by a factor of 13 compared to the full calculation. 
A similarly accurate perturbative calculation based on 
Eq. (|ll|l is still 8 times slower, which demonstrates the 
usefulness of the extrapolation to extend the calculations 
beyond their current limits. 

We now apply the extrapolation method to the ^Li 
nucleus. The numerical results have been obtained using 
the Many-Fermion Dynamics (MFD) code [l^]. Since we 
work here in a basis of slater determinants, we guarantee 
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FIG. 7: The effect of Nm on the extrapolation error. Symbols 
from right to left correspond to Nm = Nm — 6, ... , A'^^- The 
lines follow Eo = Ol^AE^ + So and Eo = a~'^AE^ + So for 
each Nm- The lines with Nm = 20 are estimated from the 
effective kinetic energies. 



a Os oscillator state of CM motion of our physical states 
by adding a CM term to equation A(Hcm -^h^) 0, 
with A = 10. This separates excited states of CM motion 
from low lying physical states. FigureElshows the results 
with Nm — 12 (i.e., 10 hfl above the lowest unperturbed 
oscillator configuration of ^Li ^^)- The basis dimension 
of the N„i — 12 calculations is 9.7 million, which is the 
largest model space published to date for ^Li. This low 
value of Nm limits what N^ values one can use in the 
extrapolation for two reasons. Firstly, the potential en- 
ergy is not negligible at small quantum numbers, so Nm 
should be at least greater than 4 (the ground state of ^Li 
has N = 2). Secondly, the accuracy of the extrapolation 
depends on the dominance of the diagonal elements of 
H, which is also weakened at small N. Thus, Eq and 
AE'^ from Nm = 4 calculations are not in line with those 
from Nm = 6, 8, 10, and 12. The extrapolation errors 
are 290 keV (CD-Bonn) and 220 keV (AV18). The CD- 
Bonn potential provides 1.1 MeV more binding than the 
AVIS potential. Thus, it is interesting to note that the 
difference between the two potentials is larger than the 
extrapolation error. 

The range of Nm that is suitable for extrapolation is 
expected to increase with Nm ■ This is confirmed in Fig.|3 
where an additional Nm = 10 calculation is made for 
comparison. It is seen that the results are much more 
linear for the case Nm = 12. Figure [T] also demonstrates 
why a higher value of Nm is likely to produce a smaller 
extrapolation error. We can estimate the lower and up- 
per bounds of the constant a, (aL,au)j in equation H12|) 
using effective kinetic energies and unperturbed ground 
state energies Eq for Nm = Nm — 2 and Nm — 6. The re- 
sults of Nm calculations are consistent with second-order 
perturbation theory, because they are roughly bounded 
by ttL and au for each Nm- The opening angle between 



the lines of slope and a^^ decreases as Nm increases. 
Hence, we expect to reduce the extrapolation error by in- 
creasing Nm- 



V. CONCLUSIONS AND OUTLOOK 

Because of the need for a large basis space to achieve 
accurate results, even light nuclei require a significant 
amount of computing resources to be investigated in the 
NCSM. We have justified and verified an extrapolation 
method for NCSM calculations. It is reliable, and can 
provide good estimates of large-space results from several 
small-space calculations. Sometimes, it may be the only 
means for getting a useful estimate of the NCSM result 
for otherwise unachievable large model spaces. 

The extrapolation formula proposed in Ref. |^ agrees 
in leading order with our Eq. (|12|l . We would like to 
emphasize that the linear scaling between Eq and AE"^ 
is based on perturbation theory. It is not an expansion 
in terms of AE'^Eq'^ , because AE^ can be much greater 
than Eq in the NCSM (see Fig.gj). The reasoning behind 
our extrapolation method is probably applicable only to 
NCSM calculations, because we have explicitly made use 
of the structure of the Hamiltonian in our derivation. The 
linear scaling between Eq and AE'^ relies on the flatten- 
ing of the diagonal elements of the effective Hamiltonian 
(dominated by the kinetic energy) as N approaches Nm- 
With this behavior, the energy denominator in equation 
(|10|l can be approximated by a constant a. 

Generally speaking, extrapolation methods depend not 
only on the structure of the Hamiltonian but also on the 
truncation scheme that is used to produce an approxi- 
mate state. Different truncation schemes may lead to dif- 
ferent scaling behavior The traditional phenomeno- 
logical shell model does not generate the structure of 
the Hamiltonian that is advantageous for our method. 
Specifically, in calculations with a core diagonal dom- 
inance is reduced since energies relative to a core are 
obtained. For realistic mean field potentials, the single 
particle spectrum does not rise as fast as an oscillator 
spectrum which itself does not rise as fast as the kinetic 
spectrum in an oscillator basis. Hence, the diagonal dom- 
inance we have in the NCSM is much stronger than the 
traditional shell model, and our method probably cannot 
be applied to the traditional shell model without modi- 
fications. On the other hand, the methods for the tradi- 
tional shell model do not necessarily apply to the NCSM 
either. In fact, it is evident from Fig. that a quadratic 
fit 9] to the Nm = 6, 8, and 10 results will yield a signif- 
icantly larger error. 

From Table |n] we have learned that for small model 
spaces there are two competing sources of inaccuracy: 
one in the perturbation theory result and the other in 
the extrapolated result. We have shown that the ex- 
trapolated result is converging to the perturbation the- 
ory result, equation (|10|) . In Table HH one sees that the 
deviation of equation H10|) and the extrapolation is al- 



8 



ready small even for relatively small values of Nm, e.g., 
Nm = 10. Since perturbation theory requires larger Nj^, 
around 16, to give estimates of similar accuracy, we can 
conclude that overall our results are probably dominated 
by errors due to perturbation theory. 

This method, like NCSM calculations themselves, is 
limited by the size of the model-space, because Ai?^ has 
to be evaluated in the full P space. Nevertheless, it has 
been shown to be a valuable tool. Its power is based on 
the much smaller dimension of the B matrix compared 
to the full matrix for the Hamiltonian operator. Calcu- 
lations for even larger model spaces will become possible 
once we make explicit use of the small dimension of the 
B matrix in our codes. 

In Figs. 01 and El we observe that the extrapolation 
error, i.e., SE = Eq — Eq, at AE^ = is quite small 
compared with differences in exact results of £q for vari- 
ous values of hfl or choices of potentials. Therefore, the 
extrapolation can also serve as a way to estimate the un- 
certainties of NCSM results arising from hil dependence 



and choices of interactions. This will be an important ap- 
plication of our method in future investigations of nuclei 
within the NCSM. 
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